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Abstract. In this note we give a simple proof of the classical fact that the 
hexagonal lattice gives the highest density circle packing among all lattices in 
M 2 . With the benefit of hindsight, we show that the problem can be restricted 
to the important class of well-rounded lattices, on which the density function 
takes a particularly simple form. Our proof emphasizes the role of well-rounded 
lattices for discrete optimization problems. 



1. Introduction 

The classical circle packing problem asks for an arrangement of nonoverlapping 
circles in E 2 so that the largest possible proportion of the space is covered by 
them. This problem has a long and fascinating history with its origins in the works 
of Albrecht Diirer and Johannes Kepler. The answer to this is now known: the 
largest proportion of the real plane, about 90.7%, is covered by the arrangement 
of circles with centers at the points of the hexagonal lattice. The first claim of a 
proof was made by Axel Thue in 1892, and then once again in 1910. It is generally 
believed however that the first complete flawless proof was produced only in 1940 
by Laszlo Fejes-Toth (see [2J, [TU] for detailed accounts and bibliography). On the 
other hand, the fact that the hexagonal lattice gives the maximal possible circle 
packing density among all lattice arrangements has been known much earlier: all 
the necessary ingredients for the first such proof were present already in the work 
of Lagrange, although he himself, while aware of the circle packing problem, may 
not have realized that he essentially had a proof for the optimal lattice packing 
in hands. In fact, the notion of a lattice has not been formally introduced until 
the work of Gauss in 1831. A detailed history and overview of these and other 
developments in the direction of the circle packing problem and its (much more 
difficult) three-dimensional analogue, the Kepler's conjecture, can be found in the 
excellent recent book of G. G. Szpiro [T2"] . 

In this note we concentrate on the lattice circle packing problem. Let us first set 
up the basic notation and describe the problem. Recall that a lattice A in M 2 is a 
free Z-module of rank two, so A = XI? for some matrix X — (xi x 2 ) S GL 2 (M), 
where the column vectors x±, x 2 of X form a basis for A and X is referred to as the 
corresponding basis matrix. The determinant of A, denoted by det(A), is defined 
to be |det(AT)|, which does not depend on the particular choice of a basis for A. 
Let us now construct a circle packing associated to A. Define the Voronoi cell of A 
to be 

V(A) = {yeK 2 :||y||<||y-:r||V:reA}, 
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Figure 1. Hexagonal lattice with Voronoi cell translates and as- 
sociated circle packing 

where we write || || for the Euclidean norm on IR 2 . In other words, V(A) is the 
closure of the set of all vectors in the real plane which are closer to than to any 
other vector of A. The area of the Voronoi cell is equal to det(A), and 

R 2 = (J V(A) + y, 

y£A 

meaning that the real plane is tiled with the translates of V(A). Moreover, as 
is clear from the definition, the interiors of these translates are disjoint. Let us 
inscribe a circle into each translate V(A) + y of this Voronoi cell by a point of the 
lattice, and write r(A) for the radius of this circle. No two such circles overlap, 
and so we have a circle packing in M 2 , called the lattice packing corresponding to 
A. The density of this circle packing is now given by 

.... area of one circle 7rr(A) 2 

A(A) = = — 

area of the Voronoi cell det(A) 

The lattice packing problem in R 2 is to maximize this density function on the 
space of all lattices. The answer has been known since the end of the nineteenth 
century (see Figure [l]) : this density function A on lattices in K 2 is maximized by 
the hexagonal lattice 

a -=G |) z2 - 

Here we will present a proof of this fact, emphasizing the particular properties of 
Ah that make it a solution to this optimization problem. 

Let us say that two lattices A and in ]R 2 are similar if there exists a real 
constant a and an 2 x 2 orthogonal real matrix U such that 

fl = aUA, 

in other words, if £1 can be obtained from A by rotation and dilation. Similarity is 
readily seen to be an equivalence relation, and it is easy to notice that the packing 
density function A is constant on each similarity class. We will prove the following 
classical result. 



REVISITING THE HEXAGONAL LATTICE 



3 



Theorem 1.1. Let A be a lattice of rank 2 m K . Then 
(1) A(A) < A (A/,) = = 0.906899 . 

wi£/i equality in |Ip if and only if A is similar to Aft. 



2. Background and standard notation 

We start by setting up some additional notation. Let B be the unit circle centered 
at the origin in M 2 . Given a lattice A, we define Minkowski successive minima 
Ai < A2 of A to be 

Xi = inf {A G M>o : A n Ai? contains i linearly independent nonzero vectors} , 

where i = 1, 2. By definition of the Voronoi cell of A, its in-radius is equal to one 
half of the distance from the origin to the nearest (with respect to Euclidean norm) 
point of A, which is precisely Ai/2, and so 

We will say that the vectors x±,X2 G A correspond to successive minima Ai,A2 if 
they are linearly independent and 

H^ill = Ai, ||x 2 || = A 2 . 

Notice that if Xi,x 2 correspond to successive minima in A, then so do ±£Ci,±a: 2 . 
From now on, when we refer to vectors corresponding to successive minima in a 
lattice in M 2 , we will always mean a pair of such vectors so that the angle 9 between 
them is in the interval [0,7r/2]. Therefore cos 9 > 0, and so 

(3) x\x 2 = Hccxlllla^ll cos ^ > 0- 

A lattice A C K 2 is called well-rounded, abbreviated WR, if its successive minima 
Ai and A2 are equal. The hexagonal lattice A^ is an example of a WR lattice with 
Ai = A2 = 1. Well-rounded lattices are very important in coding theory [1] and 
discrete optimization problems [7] ; they also come up in the context of some number 
theoretic problems, such as Minkowski's conjecture [5] and the linear Diophantine 
problem of Frobenius [5] . For a detailed study of the distribution of certain types 



of WR lattices in K 2 see [3] and [I]. In Lemma 3.7 below we show that the WR 
property is preserved under similarity, i.e. a well-rounded lattice in IR 2 can only be 
similar to another well-rounded lattice, and give a simple necessary and sufficient 
criterion for two WR lattices in R 2 to be similar. Thus Theorem |l.l| implies right 
away that only a WR lattice can maximize lattice packing density. 

Our proof of Theorem |1.1| emphasizes the importance of WR lattices. Specifi- 
cally, we first prove that A must achieve its maximum at a WR lattice, hence this 
optimization problem can be restricted to WR lattices only. Next we show that 
if A is WR, then A (A) is given by a particularly simple expression, and maximiz- 
ing it becomes an easy problem. Our argument is self-contained and requires no 
background beyond linear algebra. For further topics in the fascinating subject of 
lattice packing in dimensions two and higher see [2J, [5J, [7], [TO], and [TT]. We are 
now ready to proceed. 
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3. Properties of well-rounded lattices in R 

Our goal here is to prove that the circle packing density function on the space 
of all lattices in R 2 achieves its maximum at the hexagonal lattice. We start with 
a simple, but very useful lemma. 

Lemma 3.1. Let X\ and x 2 be nonzero vectors in K 2 so that the angle 9 between 
them satisfies < 6 < ^ . Then 

\\x 1 - £c 2 || < max{||»i||, ||x 2 ||}. 

Proof. Notice that x\x 2 > by pi. Then, since < f , 



1 ^ x\x 2 

< cos 9 - 



2 IMIM' 

and hence 

||zi-a: 2 || 2 = {x 1 -x 2 ) t {x 1 -x 2 ) = \\x 1 \\ 2 + \\x 2 \\ 2 -2x\x 2 

< ll^lP + ll^lp-II^IHI^I^maxlll^lUI^II} 2 . 

□ 

Lemma |3 . 1| readily implies that the angle between vectors corresponding to suc- 
cessive minima in a lattice cannot be < 7r/3. 

Lemma 3.2. Let AcM 2 be a lattice of full rank with successive minima Xi < X 2 , 
and let Xi,x 2 be the vectors in A corresponding to Ai,A 2 , respectively. Let 8 £ 
[0,7r/2] be the angle between X\ and x 2 . Then 

tt/3 < 6 < tt/2. 



Proof. Assume that 9 < tt/3, then Lemma 3.1 implies that 

\\xi - x 2 \\ < \\x 2 \\ = A 2 , 

which contradicts the definition of A 2 since the vectors Xi and Xi — x 2 are linearly 
independent. □ 

We can now prove that vectors corresponding to successive minima in a lattice 
in M 2 form a basis. 

Lemma 3.3. Let A be a lattice in M 2 with successive minima Ai < X 2 and let xi, x 2 
be the vectors in A corresponding to Ai,A 2 , respectively. Then Xi,x 2 form a basis 
for A. 

Proof. Let y 1 e A be a shortest vector extendable to a basis in A, and let y 2 S A be 
a shortest vector such that y 1: y 2 is a basis of A. By picking ±y 1 ,±y 2 if necessary 
we can ensure that the angle between these vectors is no greater than tt/2. Then 

0< Ill/ill < HVall, 

and for any vector z £ A with \\z\\ < \\y 2 \\ the pair y 1; z is not a basis for A. Since 
Xi,x 2 € A, there must exist integers a\, a 2 , b%, b 2 such that 



(4) (X! x 2 ) = ( Vl y 2 ) 



a x bi 
a 2 b 2 
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Let 6 X be the angle between Xi,X2, and 8 y be the angle between yi,y 2 , then 

3.2 Moreover, ir/3 < 9 y < ir/2: indeed, suppose 



7r/3 < X < tt/2 by Lemma 
9 y < 7r/3, then by Lemma 3.1 



\\Vi -1/2 II < WViW: 

however y\,y^ — y 2 is a basis for A since y 1 , y 2 is; this contradicts the choice of y 2 . 
Define 



V 



det 



K a 2 b 2/ 

then 13 is a positive integer, and taking determinants of both sides of Q , we obtain 
(5) ||xi||||cc2|| sin0 x = X>||i/i||||i/ 2 || sin^. 

Notice that by definition of successive minima, ||a;i||||a;2|| < 1 1 2/ 1 1 1 1 1 1 1 3 an d hence 
^ implies that 

v= llgjjllM sin 9 X < _2_ 2 

II2/1IIII2/2II sinS v ~ \/3 
meaning that V = 1. Combining this observation with Q, we see that 

{X1X2) (7 2 bl) =(»iWa), 

where the matrix [ 1 ^ ) has integer entries. Therefore cci,a:2 is also a basis 
\«2 02/ 

for A, completing the proof. □ 
Remark 3.1. We note that if we replace M 2 with M. d then the statement of Lemma 



3.3 is no longer true for d > 5 (see for instance [3]). 



We will call a basis for a lattice as in Lemma |3.3| a minimal basis. The goal of the 
next three lemmas is to show that the lattice packing density function A attains 
its maximum in R 2 on the set of well-rounded lattices. 

Lemma 3.4. Let A and fl be lattices of full rank in R 2 with successive minima 
Ai(A),A 2 (A) and Ai(f2), X 2 (fl) respectively. Let X\.x 2 and y 1 ,y 2 be vectors in A 
and Q, respectively, corresponding to successive minima. Suppose that X\ = y 1: 
and angles between the vectors X\,x 2 and y l5 y 2 are equal, call this common value 
0. Suppose also that 

Ai(A) = A 2 (A). 

Then 

A(A) > A(fi). 

Proof. By Lemma [3.3| x\ , x 2 and y 1 , y 2 are minimal bases for A and fl, respectively. 
Notice that 

Ai(A) = A 2 (A) = ||xi|| = ||a: 2 || 

= ||y 1 ||=A 1 (fi)<||y 2 ||=A 2 (fi). 



Then, by 

A(A) 



ttAi(A) 2 Ai(A) 2 7T 



4det(A) 4||xi||||x2||siii0 4sin( 



1 1 ~ 4\\ Vl \\ \\y 2 \\ S in6 4det(tt) 1 1 



6 



LENNY FUKSHANSKY 



□ 

The following lemma is a converse to Lemma |3.2| 



Lemma 3.5. Let A C 1R 2 be a lattice of full rank, and let X\,x 2 be a basis for A 
such that 

ll^ill = II^H, 

and the angle 9 between these vectors lies in the interval [7t/3,7t/2]. Then X\,x 2 is 
a minimal basis for A. In particular, this implies that A is WR. 

Proof. Let z E A. then z = axi + bx 2 for some a,b E Z. Then 

||z|| 2 = a 2 ||a;i|| 2 + b 2 \\x 2 \\ 2 + 2abx f 1 x 2 = (a 2 + b 2 + 2abcos9)\\x 1 \\ 2 . 

If ab > 0, then clearly ||z|| 2 > ||a:i|| 2 - Now suppose ab < 0, then again 

ll^ll^^ + ^-la&Dllzill^llzill 2 , 

since cos# < 1/2. Therefore Xi,x 2 are shortest nonzero vectors in A, hence they 
correspond to successive minima, and so form a minimal basis. Thus A is WR, and 
this completes the proof. □ 

Lemma 3.6. Let A be a lattice in K 2 with successive minima A l5 A 2 and corre- 
sponding basis vectors X\,x 2 , respectively. Then the lattice 

A W R = ^i Y 2 X ^) Z2 
is WR with successive minima equal to \\. 



Proof. By Lemma 3.2 the angle 9 between X\ and x 2 is in the interval [7r/3,7r/2], 



and clearly this is the same as the angle between the vectors X\ and j^x 2 . Then 
by Lemma|3.5| Awr is WR with successive minima equal to Ai. □ 



Now combining Lemma |3.4| with Lemma |3.6| implies that 
(7) A (A W r) > A (A) 

for any lattice AcR 2 , and ([6| readily implies that the equality in ([7| occurs if and 
only if A = A\vr, which happens if and only if A is well-rounded. Therefore the 
maximum packing density among lattices in M 2 must occur at a WR lattice, and 
so for the rest of this section we talk about WR lattices only. Next observation is 
that for any WR lattice A in K 2 , ^ implies: 

sin 9 = — — — , 
4A(A)' 

meaning that sin 9 is an invariant of A, and does not depend on the specific choice 
of the minimal basis. Since by our conventional choice of the minimal basis and 



Lemma 3.2 this angle 9 is in the interval [7r/3,7r/2], it is also an invariant of the 



lattice, and we call it the angle of A, denoted by 9(A). 

Lemma 3.7. Let A be a WR lattice in M 2 . A lattice f2 C K 2 is similar to A if and 
only ifQ, is also WR and 9(A) = 9(fl). 
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Proof. First suppose that A and f2 are similar. Let x 1,2:2 be the minimal basis for 
A. There exist a real constant a and a real orthogonal 2x2 matrix U such that 
il = all A. Let y 1 ,y 2 be a basis for f2 such that 

(l/i y 2 ) = o^Oi K2). 

Then Hy-J) = ||y 2 ||; an d the angle between y 1 and y 2 is 0(A) G [7r/3,7r/2]. By 
Lemma |3.5| it follows that y 1 ,y 2 1S a minimal basis for f2, and so f2 is WR and 
0(0) = 9(A). 

Next assume that is WR and 0(f2) = 0(A). Let A(A) and X(il) be the respective 
values of successive minima of A and f2. Let Xi, x 2 and y 1 , y 2 be the minimal bases 
for A and il, respectively. Define 

A(A) A (A) 

Zl = W) Vl ' Z2 = W) V2 - 

Then x\, x 2 and Z\, z 2 are pairs of points on the circle of radius A(A) centered at 
the origin in R 2 with equal angles between them. Therefore, there exists a 2 x 2 
real orthogonal matrix U such that 

and so A and are similar lattices. This completes the proof. □ 

We are now ready to prove the main result. 

Proof of Theorem \1.1\ The density inequality ^ says that the largest lattice pack- 
ing density in W 2 is achieved by some WR lattice A, and (JsJ> implies that 

(8) A(A) 



4 sin 9(A)' 



meaning that a smaller sin 9(A) corresponds to a larger A(A). Lemma 3.2 implies 
that 9(A) > 7r/3, meaning that sin 9(A) > V3/2. Notice that if A is the hexagonal 
lattice 

A * = i) z! ' 

then sin 9(A) — v3/2, meaning that the angle betw een the basis vectors (1,0) 
and (1/2, y/3/2) is 9 = tt/3, and so by Lemma 



3.5 



this is a minimal basis and 
9(A) — tt/3. Hence the largest lattice packing density in M 2 is achieved by the 
hexagonal lattice. This value now follows from (Js]> . 

Now suppose that for some lattice A, A (A) = A(A^), then by Q and a short 
argument after it A must be WR, and so 

A(A) = — ^— = A(A h ) 



4 sin 6»(A) v 4sin7r/3 
Then 9(A) = tt/3, and so A is similar to by Lemma |3.7| This completes the 



proof. □ 
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